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We explore lepton-flavored electroweak baryogenesis, driven by CP-violation in leptonic Yukawa 
sector, using the 7 — u system in the two Higgs doublet model as an example. This setup generically 
yields, together with the flavor-changing decay h —> Tu, a tree-level Jarlskog-invariant that can drive 
dynamical generation of baryon asymmetry during a first-order electroweak phase transition and 
results in CP-violating effect in the decay h — rr. We find that the observed baryon asymmetry 
can be generated in parameter space compatible with current experimental results for the decays 
h > Tu, h > TT and T > uy, as well as the present bound on the electric dipole moment of the 
electron. The baryon asymmetry generated is intrinsically correlated with the CP-violating decay 
h — 77 and the flavor-changing decay h —> Tpu, which thus may serve as “smoking guns” to test 


lepton-flavored electroweak baryogenesis. 


Introduction. Explaining the origin of the baryon 
asymmetry of the universe (BAU) is a forefront challenge 
for fundamental physics. The BAU is characterized by 
the baryon density ng to entropy s ratio 

NB 


=e S (8.61 + 0.09) x 107" [I]. (1) 


According to Sakharov [2], generation of a non-vanishing 
Yp requires three ingredients in the particle physics of 
the early universe: non-conservation of baryon number 
(B); C- and CP-violation (CPV); and out of equilibrium 
dynamics (assuming CPT conservation). While the Stan- 
dard Model (SM) of particle physics contains the first 
ingredient in the guise of electroweak sphalerons, it fails 
with regard to the remaining two. Physics beyond the 
SM is, thus, essential for successful baryogenesis. 

Electroweak baryogenesis (EWBG)[3] is among the 
most theoretically well-motivated and experimentally 
testable scenarios, as it ties BAU generation to elec- 
troweak symmetry-breaking (see [4] for a recent review). 
Extending the SM scalar sector can lead to a first order 
electroweak phase transition(EWPT), thereby satisfying 
the out-of-equilibrium condition. Addressing the second 
Sakharov criterion requires new sources of CPV, as the 
effect of CPV in the SM Yukawa sector is suppressed 
by the small magnitude of the Jarlskog invariant associ- 
ated with the Cabbibo-Kobayashi-Maskawa (CKM) ma- 
trix and by the small quark mass differences relative to 
the electroweak temperature, Tew ~ 100 GeV. 

It is possible that an extended Yukawa sector may rem- 
edy this SM shortcoming. A particularly interesting yet 
unexplored possibility involves the leptonic Yukawa in- 
teractions. Phenomenologically, the report by the CMS 
collaboration of a signal for the charged lepton flavor 
violating (CLFV) Higgs boson decay h > Ty (2.4 o 


significance) [8] hints at a possible richer leptonic Yukawa 
sector [9], though the ATLAS collaboration observes no 
evidence for this decay mode[I0]. Should an extended 
leptonic Yukawa sector exist, then the accompanying new 
CPV phases may provide sources for EWBG that do not 
suffer from the suppression associated with SM quark 
Yukawa sector. 

Motivated by these considerations, we study the via- 
bility of “lepton flavored EWBG”, a scenario that relies 
on both CLFV and leptonic CPV. For concreteness, we 
use a variant of the type III two Higgs doublet model 
(2HDM) [I5] with generic leptonic Yukawa textures [I6] 
and focus on the T — u families as an example. For a 
representative choice of Yukawa texture, we derive the 
CPV source for the EWBG quantum transport equa- 
tions [II] [13] in terms of the relevant Jarlskog invari- 
ant, ImJ,4. We then solve these equations, which encode 
the dynamics of CLFV scattering during the electroweak 
phase transition, and obtain the BAU as a function of the 
Yukawa matrix parameters. We also show that the same 
ImJ,4 also generates a CPV coupling of the Higgs boson 
to T leptons at T = 0, parameterized by a CPV phase 
¢,. Measurements of CPV asymmetries in h > ttT~, 
as discussed in Ref. [14], would provide a test of this 
baryogenesis mechanism. Taking into account present 
constraints from measurements of T (h > 7+77) and lim- 
its on T(t > py) we find that a O(10°) determination of 
r would probe this scenario at a significant level. 
Model Setup. Since our focus is on CPV in the u — 
T sector which has no mixing with the first-generation 
charged leptons, we assume the scalar potential to be CP- 
conserving with parameters chosen to generate a strongly 
first order EWPT [25]. The particle spectrum consists of 
two CP-even neutral scalars (h, H), the neutral CP-odd 
A, and a pair of charged scalars H*, with the lighter h 


taken to be SM-like. The SU(2), x U(1)y invariant weak 
eigenbasis lepton Yukawa interaction is 


a = -L [Yi i781 + Yo i; B2] eh + h.c., (2) 
where ©; 2 are the two Higgs doublets with the same 
hypercharge, Lê and eh are left-handed lepton doublet 
and right-handed lepton singlet in weak basis, with the 
family index i,7 = 2,3. Then we can uniquely define a 
Jarlskog invariant as the imaginary part of [I9] [20]: 
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with the power of Yukawa coupling (or mass parameter 
of fermions) product being two. Here va = V2(®°) is 
vacuum expectation value (vev) of neutral Higgs fields, 
Hab is the coefficient of Pİ ®, in the potential, and the 
trace is taken over flavor space. J,4 is normalized to be a 


dimensionless quantity by dividing a factor v?u#?, where 


1 : 
iy = 5 (22 — M11) sin 28 + p12 cos 28 (4) 
is a quardratic Higgs coupling defined in “Higgs ba- 
sis” (15) (20): Ay = cos B®, + sin Bo; Hə = — sin B®, + 
cos B®; (H?) = v/V2 = 174 GeV; and (H9) = 0. 
The mass matrix for fermions is defined as 


M = (vi Yı + v2Y2)/ V2 (5) 


in the weak basis, with a determinant of MİM or M 
close to zero (since m, % 0). For illustration, we choose 
a texture with Yj22 = Yj23 = 0, with j = 1,2. This 
immediately yields 


Im (Ja) = —Im (¥1,32Y5'32 + Y1 ,33Y2 33) (6) 
or 
Im (Ja) = —Im (¥1,32¥5'39) , (7) 


with a further assumption Y1 33 = Yo.33. The diagonal- 
ization condition |M32|?+|M33|? = m2 immediately gives 
|M32| < m+, and fixes the value of |Y1 33| = |Y2,33|. Since 
the proposed mass texture is not invariant under basis 
transformation of ®; and ®2, tan 3 = v2/vı becomes an 
independent parameter (similar to what happens in type 
II 2HDM). Thus this setup contains five relevant and in- 
dependent parameters: tan 8, a (the mixing angle in the 
CP-even Higgs sector), |Y2,32], r32 = [Yi 32|/|Y2,32] and 
Im (Ja). 

In the mass basis for both fermions and Higgs bosons, 
the r Yukawa interaction is then parameterized as 


1 
— g eTRlh(mrsp—a + N;7Cg—a) 
+H(m-cg—a — NrrSp-a) +iAN-7| + h.c., (8) 


where 6 — a is invariant under the basis transformation 
in Higgs family space [I8]. The SM-like Higgs boson h 


receives two contributions to its coupling. The first one 
results from its H? component which is aligned with the 7 
mass. Another one is related to its HÌ component which 
is proportional to N,,, the Yukawa coupling of H9 with 
T leptons, with 


v? uty Re(Ja) — 2u m? 
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Re(N,-) = 
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The CLFV interactions are completely controlled by the 
Yukawa coupling of H9, N, 


N+ : 
— —, TEER (Cp—ah — S6-aHĦ + iA) +h.c., (10) 
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With tan 8 = 1, the expression in terms of weak basis 
parameters is given by 


M33 


N, = et? ede 
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Here ô is an un-physical phase undetermined in the diago- 
nalization procedure which can be removed by field redef- 
inition. For later convenience, we also have for tan 8 = 1 


Re (Ja) = 


(\Y2,321? — [Y1 ,321°) + 2B: (2) 
1 


NI = 
N 


Finally the charged Higgs Yukawa interactions are gov- 
erned by —V2/vHt4 Nijeh +h.c.. 

Given the four free parameters left for describing tree- 
level Yukawa interactions of the u— rT system, we present 
various phenomenological results (e.g., h => TT, Tu and 
T — uy constraints) and the BAU analysis in terms of 
the effective h?r coupling [87] (see Fig. 
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7 [Re(y-)T7 + Im(yr)TiysT]h (13) 
with benchmark values assigned to r32 and 6 — a. Here 


Re (yr) = 58-a + "Re (N;r) ; 


e—a Tm (Np-) . (14) 
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Then, the condition |M32| < m, imposes a constraint 
at the (Re(y-), Im(y-)) plane, allowing only a circular 
region centered at (Re(y,) = 8g-a + ¢g—a(1 + r3,)/(1— 
73>), Im(y,) = 0) with a radius 2|cg_ar32/(1 — r3,)|. At 
its boundary, we have M33 = 0 and hence N,, = 0. 
For r32 = 1, Nr- is purely imaginary, yielding a vertical 
line at Re(y-) = sg—a. In Fig. |1| we present results in 
two representative cases: r32 = 0.9 and r32 = 1.1, with 
B=a— > = 005. 


h — TT constraints. The decay width for h > TT is 
given by 


ee V2Grmy,m2 
pt = ly? 
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Experimentally, the ATLAS signal strength is wAtpas = 
LAs eee while CMS favors a smaller one uGhiygs = 
0.78 + 0.27 [35]. We take a x? analysis at 95% C.L. for 
these two measurements, assuming a Gaussian distribu- 
tion for both and neglecting their correlations. Appar- 
ently, the allowed parameter region should be a circular 
band at the (Re(y-), Im(y-)) plane, as is indicated by two 
green dashed curves in Fig.{1} A future determination of 
this coupling that agrees with the SM value within +10% 
is plotted as a curved blue band. 

h — Tp constraints. The lepton flavor-changing decay 
width is given by 
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[Nerul = (16) 
Theoretically, a sizable Br(h — Ty) requires a small 
|M32| (see Eq. (11)). ATLAS sets an upper limit on its 
branching ratio, Br(h + Ty) < 1.43%, at 95% C.L. [10], 
while CMS gives a best fit Br(h > Ty) = 0.847339% 
as well as an upper limit Br(h —> Tu) < 1.51% at 
95% C.L. [8]. In Fig. the current ATLAS upper 
limit 1.43% as well as two projected ones, say, 1% and 
0.5%, are shown as dashed curves, in the two cases with 
r32 = 0.9 and 1.1. The circular boundaries of the brown 
regions correspond to vanishing M33 or N-,,, yielding 
Br(h > Tu) = 0. 

T — py constraints. Non-vanishing N;, may also 
contribute to the rare decay T + py, via one-loop neutral 
and charged Higgs mediated diagrams and two-loop Barr- 
Zee type diagrams [32] [83]. Explicitly, one has 


T,aG?.m> 
3274 


where T, = (290.3+0.5) x 10~1°s [30] is the 7 lifetime and 
C7L/r are the Wilson coefficients of the dipole operators 


Br(t > uy) = (ICrx|? +|CFpl), (17) 


= — m, ño” (PrP yes (18) 
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in the Hamiltonian —Gp[C7,Q¥ + CrrQ?|/V2B]]. In 
our setup, C7z and Cyr are proportional to N7,, and 
Nur, respectively, yielding a vanishing Cyr. The cur- 
rent experimental limit is Br(r > py) < 4.4 x 1078 
(90% C.L.) [29]. The allowed parameter regions are de- 
noted in gray in Fig.|1} Obviously there exists a positive 
correlation between the experimental constraints from 
Br(h + Tu) and Br(v > py), though the relevant new 
physics contributions result from tree- and loop-levels, 
respectively. 

Electric dipole moments. Null results from experi- 
mental searches for the electric dipole moments (EDMs) 
of the neutron, neutral atoms, and molecules in gen- 
eral place stringent limits on new sources of CPV. In 
the present instance, the electron EDM (de) provides 
the most significant probe of Im (Ja) or Imy,, given 
the bound obtained by the ACME collaboration using 


ThO/[42]. In our setup, the dominant contribution to elec- 
tron EDM results from h—mediated Barr-Zee diagram 
with a 7 lepton loop, because of non-vanishing Im y+. We 
find |d./e| ~ 1.66 x 107?’ |Im y,|cm, yielding a bound of 
|Imy,;| < 5.2. As indicated in Fig. |1| this bound is an 
order of magnitude larger than what is required to ac- 
count for the observed BAU (see below). We also note 
in passing that CPV in the scalar potential, associated 
with a different Jarlskog invariant, will lead to mixing 
between the CP-even and CP-odd scalars. The result- 
ing EDM contributions can be considerably larger (for 
a given value of the relevant Jarlskog invariant). For a 
recent analysis, see Ref. |36). 


Electroweak baryogenesis. The first order EWPT 
proceeds via bubble nucleation. CPV scattering from the 
bubble walls generates a net left-handed fermion density 
nz, whose diffusion ahead of the advancing wall biases 
the electroweak sphalerons into producing a net baryon 
number density, ng. The expanding bubbles capture and 
preserve this density if the sphaleron processes are suf- 
ficiently quenched inside the bubbles. We compute nz 
from a set of quantum transport equations, derived from 
the equations of motion for Wightman functions arising 
in the closed time path formulation of non-equilibrium 
quantum field theory by expanding in gradients of the 
bubble wall profile and chemical potentials (see Ref. [I3] 
for pedagogical discussions). As with earlier work, we will 
employ the “vev insertion approximation”, which pro- 
vides a reasonable estimation of the CPV sources (see 
Ref. for a discussion of theoretical issues associated 
with the computation of these sources), and work in 
fermion weak basis (i.e. approximate fermion mass basis 
during EWPT) in this section. Since the weak sphaleron 
rate [ys [28] is much smaller than the rates for diffu- 
sion and particle number changing reactions that govern 
nz |40|, we first solve for this density and substitute the 
result into the equation for npg. 


For simplicity, we neglect bubble wall curvature [27], 
so that the quantities entering the quantum transport 
equations depend only on the coordinate in the bubble 
wall rest frame Z = z + vyt with v, being the wall ve- 
locity, Z > 0 corresponding to broken phase and Z < 0 
for unbroken phase. Since non-zero densities for the first 
and second generation quarks as well as for the bottom 
quark are generated only by strong sphaleron processes, 
the following relations hold: Qı = Q2 = —2U = —2D 
—2C = —2S = —2B, where Q,; denotes the density of 
left-handed quarks of generation k and U, D, etc. de- 
note the corresponding right-handed quark densities. In 
addition, Lı = Lə = er % 0 since the corresponding 
leptonic Yukawa interactions are negligible compared to 
those retained in our choice of Yukawa texture. Local 
baryon number density is also approximately conserved 
on the time scales relevant to the reactions that govern 
nz, so that 52 (Qi+Ui+ Di) = 0. The resulting trans- 
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FIG. 1. Theoretical and phenomenological constraints on the Higgs-r Yukawa couplings in Eq. (13). The inner parts of circular 
regions satisfy the diagonalization constraint |M32| < m+ for two representative choices of r32, with the outer boundaries giving 
vanishing Br(r > wy) and [(h > Ty). The r32 = 0.9 and r32 = 1.1 regions are separated by the vertical dashed line at 
Re(y-) = sin(0.05 + $) ~ 1. Brown regions correspond to non-vanishing T (h > ty), with different representative values (1%, 
0.5% and 0%) indicated by circular dashed lines. For r32 = 1.1, the ATLAS 95% C.L. upper bound of 1.43% is shown, while 
for r32 = 0.9 a maximum BR of 1.41% can be achieved within the theoretically allowed region. Upper limits on [(h > tr) 
(95% C.L.) and Br(r > py)(90% CL) are given by the green and grey regions, respectively. The region inside the green dashed 
lines gives the Higgs signal strength 4177 allowed region at 95 % C.L. without assuming a specific Yukawa texture. The inner 
light-blue band labelled |y-| = 1 + 0.1 corresponds to the region with a more SM-like h?r coupling. The region giving the 
observed BAU is indicated by the horizontal pink bands assuming |AG| < 0.4) for 6—a— 5 = 0.05 as discussed in the text. The 
other relevant parameters are fixed to be my = ma = my+ = 500GeV, vw = 0.05, Lw = 2/T, Da = 6/T and T = 100GeV. 
To guide the eye, the argument of y- is indicated with red-dotted lines. Note, the calculation of baryon asymmetry outside the 
circular regions could be unreliable due to the breaking of perturbative “mass insertion”. 


across the bubble wall. Furthermore Iss ~% 16a47 is 


the strong sphaleron rate |24|; Tm: is the two body top 


port equations are 


H— \ 
OQ = Tm (Er — Qs) +Te(Sr — EH — SQ) relaxation rate [I3]; and T;/+ is the t/r Yukawa induced 
+20 ss0ss, three body rate [26]. After solving for the densities in 
3 LH =Ti(Er — £y — £93) +D- (Eng — Err — ÊH) Eqs. (19), we obtain ng = >>,(Qi; + Li) 27] and ng, 
which is a constant in the broken phase: 
—20 én, 
3 L5 = -Imr (EL — Err) — Ur (Ets — Erg — Ên) np= 3D ws | ny(Z)e-7dz , (21) 
gcoPV DgA+ 
Outer = T, (ên + Ern — Era) + Dmr (êr — Ern), where Ty, œ~ 120a T PBS] and à} = (vw + 


v2, + 15r'wsD4)/(2D4). 
Assuming a fast Tr diffusion [41], we solve the trans- 
port equations perturbatively at the leading order of Ty a 
(19) jt, Py* andTy,'. We have further neglected Pmr in the 
final result as it is generally small compared with Imz; 
then ng is proportional to Im(y,-) with no dependence 
on Re(y,). One important remaining parametric uncer- 
tainty is the difference of 6(Z) in the broken and symmet- 
ric phases (= A£) since the CPV source term and thus 
np are both directly proportional to it. Here we take its 


ô T! T —Pme(Er Qs) 
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where ds, = Er + 9EB — 2€Q3, a = Na/ka, with ka be- 
ing the statistical weight [I3] associated with the number 
density nq of species “a” and 0, © Uwe — Dees with 
Da being the diffusion constant [22] from the diffusion 
approximation. The CPV source terms are 


a dß(z 
sory = gorv CEE Ima) 5 Og 
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where Z is a momentum-space integral that depends on 
the leptonic thermal masses (see Ref. [23]) and d8/dz 
characterizes the local variation of tan 6(Z) as one moves 


maximum magnitude to be 0.4 and vary it to obtain the 
bands in Fig. |1| where the upper and lower bands give 
opposite signs of BAU resulting from the unknown sign 
of AS. Imposing the condition |M32| < m, as discussed 
above then restricts Re(y,) to the region of overlap be- 
tween the pink bands and the two circular regions. 


Results and collider probes. Combining the analy- 
ses above, we find that there exist parameter regions in 
Fig. [1] where the observed BAU can be explained without 
violating current experimental bounds. These regions 
are characterized by |Im (y+)| 2 O(0.1), corresponding 
to |Im(J4)| 2 O(10~°), or |¢,| > O(10°). As indicated 
above, the present EDM upper bounds on these CPV pa- 
rameters are roughly an order of magnitude larger than 
the BAU requirements. The next generation searches for 
neutron, atomic, and molecular EDMs that plan for order 
of magnitude or better improvements in sensitivities may, 
thus, begin to probe the BAU-viable parameter space. 

Alternatively, collider measurements of the CP prop- 
erties of the h77 coupling may also test this scenario. 
For example, a recent study shows that use of the p- 
meson decay plane method or impact parameter method 
at the LHC may allow a determination of ¢, with an 
uncertainty of 15°(9°) with an integrated luminosity of 
150fb~'(500fb~'), or ~ 4° with 3ab~' [B7]. At Higgs fac- 
tories, 6, could be measured with an accuracy ~ 4.4°, 
with a 250GeV run and lab * luminosity [88]. There- 
fore, the collider measurements of the CP-properties 
of the htr coupling complement the measurements of 
h > Tu or T > py, which constrain more the parameter 
regions with relatively small |Im y,|, or |¢,|. 
Discussion. We stress that N+, or Nr are not involved 
in any Jarlskog invariants here. Indeed, for the Yukawa 
texture considered in this study, we have Ny, = Nur = 
0. Consequently, there is only one independent phase in 
the 7 — system (apart from an overall common phase), 
with the corresponding Jarlskog invariant defined by the 
imaginary part of Eq. {3}. This implies that CP-violation 
in this setup does not affect Yru Or Yur- 

With a more generalized texture, the situation can be 
different [20]. For example, one can define a new Jarlskog 
invariant [20] 


Im (Jp) ~ Im (APP Ny NrpNyurmr) (22) 


with \#® being the coefficient of the term Hi HoH! Ho in 
Higgs potential. Here the value of N,,,, need not be small, 
since Hə does not directly contribute to SM fermion mass 
generation. Nonetheless, the effects of Im (Jg) on both 
T = 0 phenomenology and CPV dynamics at finite T 
arise a loop level, since four instead of two Yukawa cou- 
plings (or fermion mass parameters) are involved. Ex- 
plicitly, Im(Jg) can be probed via the interference be- 
tween the tree-level and one-loop diagrams in the de- 
cay of h > pr. Similarly, Im (Jg) can contribute the 
BAU generation at one-loop level, as it enters the one- 
loop, self-energy correction of fermions that contribute 
to spacetime-dependent vev scattering, yielding a loop- 
level CPV source. For similar reasons, its contribution 
to muon EDM typically requires extra mass insertions, 
hence being negligibly small. A full exploration in this 
regard will be deferred to future work. 


[Note added] While this article was being finalized, the 
paper appeared, which partially overlaps with this 
one in discussing the correlation between electroweak 
baryogenesis and the 7 — u Yukawa structure in the 
2HDM. However, in contrast to Ref. [45], we having ex- 
plicitly show how the baryon asymmetry generated in- 
trinsically correlates with the CP-violating decay h —> 
TT, using a Jarlskog invariant uniquely defined at tree 
level. As for the complex phases of Higgs couplings with 
7 and u leptons, which were addressed in Ref. [45], we 
note that they contribute only to Jarlskog invariants de- 
fined as a product of at least four Yukawa couplings (or 
fermion mass parameters), and hence play a sub-leading 
role in this exploration generically. 
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